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Abstract: In this work we present an extended version of the Friedrichs Model, 
which includes fermion-boson couplings. The set of fermion bound states is cou- 
pled to a boson field with discrete and continuous components. As a result of the 
coupling some of the fermion states may become resonant states. This feature 
suggests the existence of a formal link between the occurrence of Gamow Reso- 
nant States in the boson sector, as predicted by the standard Friedrichs Model, 
with similar effects in the set of solutions of the fermion central potential (Gamow 
fermion resonances). The structure of the solutions of the model is discussed by 
using different approximations to the model space. Realistic couplings constants 
are used to calculate fermion resonances in a heavy mass nucleus. 

PACS: 03.65 Nk, 03.70 +k 



1 Introduction 

The conventional approach to the quantum mechanical description of the 
nuclear many body problem is, of course, the shell model pQ. The task 
of constructing many-body nuclear wave functions, starting from the shell 
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model treatment of two-body interactions, is formidable because of the ex- 
tremely large dimension of the basis, and only very recently some impressive 
steps towards that goal have been achieved [2]. Alternatively, one may think 
of physical oriented approximations to the problem, based on the use of effec- 
tive one-particle and collective degrees of freedom and their couplings [3J, H] . 
Central to this picture is the role play by one particle fermion bound states, 
collective boson excitations (either vibrational or rotational) and particle- 
boson couplings [I]. 

In both approaches, either the shell model or the one-particle plus col- 
lective excitations, single-particle states are obtained as solutions of nuclear 
central potentials, like the nuclear harmonic oscillator or the Woods-Saxon 
potential [3]. From the mathematical point of view these are solutions of 
the Schroedinger equation with a central potential with volume, surface and 
spin-orbit interactions [5] . Since the nucleus is a system with a finite density 
of one-particle states, one is interested in one-particle bound state solutions. 
Shell model basis are thus composed by these single-particle states j6]. Since 
the one-particle potentials are of finite range, single-particle resonant states 
are obtained by imposing purely outgoing boundary conditions [TJ. 

Tore Bergreen [8j [11] was one of the first physicist to show that a single- 
particle basis may contain not only bound states but also some single-particle 
resonances. Lately, the work of many others nuclear structure physicist fol- 
lowed [HJCE21E32] and by now the subject is rather well established and several 
nuclear structure calculations have been performed in Bergreen's basis to de- 
scribe nuclear structure properties in the continuum, decay widths of particle 
states and vibrational states, nuclear clusters and their decays [T9] . 

Since the calculations involved a rather large number of states, even in 
the simplest versions, one may be tempted to study the effects of couplings 
between bound and resonant states in the context of solvable models, in order 
to extract signatures of possible dominant mechanisms. 

In the last 40 years, the use of solvable models has contributed signifi- 
cantly to the understanding of the nuclear dynamics, in spite of its complex- 
ity. Although solvable models are popular in the context of bound single- 
particle and collective states [EH Ej, their use to describe collective and 
single-particle properties in the continuum is less frequently found in the 
literature. 

The Friedrichs model [20] is a suitable mathematical model to describe 
the coupling with the continuum. The model has the advantage of being 
exactly solvable and it can be applied to various systems of physical interest 
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[2TI 1221 [221 123]. In its simplest version the model deals with the coupling 
between discrete and continuous components of a boson field |20j. More 
sophisticated versions of the Friedrichs model have been presented in [2T| 
[221 12H]. It is the only completely solvable model which describes resonance 
phenomena [27J. A recent presentation of the model, in connection with the 
study of resonances in field theory, can be found in [23]- In this work we 
aim at the extension of the model, in order to include the coupling between 
fermion and boson fields. 

The advantage of using an exactly solvable model, like the Friedrichs 
model, to study the coupling with the continuum is rather obvious, since 
until now the description of these couplings are limited to purely numerical 
approaches. In the standard approach, like the one advocated in the litera- 
ture [9j H31 [13] , the occurrence of resonances in the fermion sector is a result 
of the boundary conditions imposed to the diagonalization of finite range 
central potentials. Purely outgoing boundary conditions lead to Gamow res- 
onances; some of these states, depending upon the imaginary part of their 
energies, may be included in single-particle basis, the so-called Berggren ba- 
sis, together with fermion bound-states. In this hibrid basis one may attempt 
to include residual two-body interactions. As shown [13] in the energies of 
some particle-hole excitations may become complex. The same effect shows 
up at the level of attractive two-particle configurations [211 [17] , where reso- 
nances may appear as a consequence of the coupling between configurations 
of fermion pairs in bound and resonant states. Finally, it is worth mentioning 
that the observed properties of atomic nuclei in the drip line may be strongly 
linked to the coupling between bound and resonant states. Particularly, we 
refer to the broadening of the single particle energies and to the appearance 
of a decay width in single-particle orbits, which may be strongly influenced 
by the coupling to continuous particle-hole excitations [T7J [T5] . A possible 
mechanism to explain for this class of phenomena is the coupling between 
bound fermion states and discrete and continuous bosonic degrees of freedom. 

Since the standard Friedrichs model describes the coupling between the 
discrete and the continuous components of a single boson field, it seems 
natural to include fermions and their couplings with the boson field and take 
advantage of the solvable structure of the model. In this respect, the results 
of the extended Friedrichs model may be used to approximate some specific 
features of fermion states, like their decay width, in situations where the 
coupling between fermions and bosons indicates the existence of effects due 
to the continuum |18j . 
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The present paper is organized as follows: In Section 2, we present the 
needed mathematical details about the model. In Section 3, we explore the 
structure of the solution obtained under different approximations. We have 
applied a compact formalism to solve the system of coupled equations yielded 
by the model. The results of applications, with reference to a realistic case, 
are presented and discussed in Section 3. Conclusions are drawn in Section 
4. 

2 The Model 

For the benefit of the reader, we shall begin with a review of the basic ideas 
of the Friedrichs Model, which is an exactly solvable model where a discrete 
boson state interacts with a continuum spectrum of boson excitations. The 
situation resembles very much that of the nuclear excitations of even mass 
nuclei near the drip line, where few bound states may strongly interact with 
low-lying continuum states. This model provides the insight about the mech- 
anism leading to the appearance of resonances in the boson sector, and it has 
the clear benefit of being solvable. 

Next, we shall review briefly the properties of resonant solutions of the 
Friedrichs model and then we shall proceed with the description of the exten- 
sion which we propose in this paper. It consists of the inclusion of fermion 
bound states which interact with the boson resonant state. Since the ex- 
tended model is also solvable, one might also get some insight about the 
effects of the coupling of boson and fermion states in the vicinity of the 
nucleon drip line. From the mathematical point of view the extension of 
the Friedrichs model amounts to the simultaneous treatment of fermion and 
boson sectors of a Hamiltonian which includes resonant states in its spec- 
trum. From the point of view of physics it represents the solvable limit of 
nucleon-boson interactions in Berggreen's basis [T8] . 

In spite of some formal steps of the derivation, the present discussion is 
motivated by the need to build a schematic, albeit solvable, model which 
may be use to approach dominant features of the particle- vibration coupling 
in the continuum. 

Further details about the formalism concerning the standard Friedrichs 
model can be found in [29J EI] and references quoted therein. 



4 



2.1 The Standard Friedrichs Model 



The simplest form of the Friedrichs model [20J includes a free Hamiltonian Ho 
with a real positive continuous spectrum (u > 0) and a discrete eigenvalue 
(a>o > 0) imbedded in the continuous spectrum. An interaction is acting 
between the continuous and discrete parts of H Q by means of a potential 
V. As a result of the action of V, the bound state of H is dissolved in 
the continuous and a resonance is produced. The spectrum of the total 
Hamiltonian H = Ho + AV^, where A is a real coupling constant, is purely 
continuous and coincides with the real positive semiaxis. We write for the 
free Hamiltonian, H , the expression 

POO 

H = u \l)(l\+ u\u)(u\du. (1) 
Jo 

The interaction, V, is given by 

poo 

V= [V*(u)\u)(l\ + V(u)\l)(u\]dw, (2) 
Jo 

where, V(lo) is called the form factor [22, 30j. 

The most general form of the state vector tp, in the basis spanned by 
{|1), \u)}, can be written as 

POO 

ip = a\l) + / ip(uj) \uj) du. (3) 
Jo 

Note that the set {|1), \u)}, u G [0, oo), of generalized eigenvectors of H , is 
a complete set. Thus, each vector state i)j has a component on the bound 
state |1) of H and an infinite number of components, in the continuous part 
of the spectrum of Ho. The action of V on ip is given by 

(POO \ POO 

J duV(u)ip(u)J |1) + a J duV*(u)\u). (4) 

To obtain the explicit form of that is to determine the amplitudes a 
and (p(u), we shall use the technique presented in |3lf*l In the notation of 
[31], the reduced resolvent, rj(z) is given by (see also [201 [30]) 

V(z) = (llj^-Jl) . (5) 



1 See also EH, section 3.5. 
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By making adequate assumptions on the form-factor V(u), rj(z) has the 
following properties [2%] : 

i. ) The function rj(z) has no singularities in the complex plane other 
than a branch cut coinciding with the positive semi-axis, i.e., the continuous 
spectrum of H. 

ii. ) The function rj(z) can be analytically continued through the cut. 
These are the functions rj±(z) 

W = 1 77 Tr^ 1= ^o-^-/ — — • 6 

H — uj ± zO Vo z — ± zO 

These extensions have poles located at the points z = E R — i^, for r]^(z), 
and Zq, for 77+ (V), with > and T > 0. 

Although resonance poles may, in principle, have arbitrary multiplicity 
[32] . depending on V(u), we shall assume hereafter that the resonance poles 
are simple. 



2.2 Resonance or Gamow states. 

Poles of (jSJ) appear into complex conjugate pairs, each of them represents a 
resonance. Then, Gamow vectors are solutions of the eigenvalue equation, 

(H - x) V(x) = . (7) 

with the total Hamiltonian H, whose eigenvalues coincide with resonance 
poles [3T] . 

Since |1) and \u) form a complete system, we must have: 

POO 

= a(x) |1) + / ip(x,u) \u) did. (8) 
Jo 

If we apply H to (jHJ), we obtain the following system of equations: 

POO 

(u - x)at(u) + A / ^(x,^)y*(a;)(ic<; = 0, (9) 

(u-i)^(i,w) + Ay(w)a(w) = 0. (10) 

To solve this system, we write a(u) in terms of if)(x,u) using (|T0|) and carry 
the result to ([9]). We obtain an integral equation which gives one solution of 
the form: 
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M, + ( x ) = \ x ) + \V*(x)^ r - \\l) + \ V{U) \u)<kA . (11) 

The bracket with an arbitrary vector of the form if), gives a function 

which is analytically continuable from above to below through the spectrum 
of the Hamiltonian H . This continuation has a simple pole at zq so that we 
can write on a neighborhood of zq\ 

^ + {z) = -^ + o{z). (12) 

Z — Zq 

From ([7]) and (TT2l) . we get 

(H - z)m + {z) = — (H - z)C +(H - z) o(z) = , (13) 

Z — Zq 

which gives 

(H - zq) C = =► EC = z C . (14) 

This shows that the residue C of *&+{z) at the pole z$ coincides, save for an 
irrelevant constant, with the decaying Gamow vector |/o)o To calculate its 
explicit form note that (1121 . on a neighborhood of Zq, has the form: 

. , constant f. . , f°° V(u) . , 1 , , 

« < 1) + A / — — \u) du> + RT , 15 

where RT stand for "regular terms" . Now, let us use the Taylor theorem to 
write: 

1 1 Z — Zn , . . . 

+ o{z). (16) 



z — u) + iO z — uj + iO (z — uj + zO) 
By replacing fflBT) in ([15"]) , we get 

. , constant f . . , Z" 00 Vfu;) , , . 



2 Each resonance is characterized by a pair of complex conjugate poles, zq = Er — iT/2 
and Zq = Er + iT/2, in the continuation of rj(z) from above to below, T}-{z), and from 
below to above, r/ + (z), respectively. Both are eigenvalues of H and the corresponding 
eigenvectors are the decaying Gamow vector, |/q), with Zq, and the growing Gamow vector, 
/o), with z%. 
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Therefore, up to an irrelevant constant, we conclude that 



The system given by equations ^ and ( JTOl has another solution that can 
be analytically continued in the upper half plane.^This solution gives, using 
the same technique, the growing Gamow vector \f ): 

i/.)Hi)+ r . Ay ' M n M^- 

Jo z -u-iO 

So far, we have introduced the simplest and more standard version of the 
Friedrichs model in a fashion accessible to the reader. This will facilitate the 
comprehension of the more complicated version under consideration in the 
next subsection. 



2.3 The extended Friedrichs Model. 

The present version of the Friedrichs model, is an extension of the previously 
introduced standard Friedrichs model of Section 1.2 and it includes : 

i) The unperturbed fermion and boson Hamiltonian Hi 

POO 

H I = u \l)(l\+ duju\uj)(u\ + S^c k \k)(k\ , (20) 

where the index k runs out the set of Fermion kets | k) . 

ii) The interaction between fermions, \k), and the discrete boson, |1), 
Hamiltonian Hif. 

#/7 = X>dM><*l + ^IO<MI] • (21) 

k,l 

hi) The interaction between fermions and the boson field, {|u;)}, Hm: 

POO 

Hiii = 5^ / du [f kJt ( U )\k,u)(l\ + WM^M] ■ (22) 

The standard Friedrichs model includes the boson-boson coupling V de- 
scribed by equation (T5]). This coupling V can be generalized to include 
fermion-boson interactions in the following manner [U [33] 
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POO 

H IV = Y J / du;[g kk ,(u J )\k,l)(k',u J \+ g* kk ,(u)\k',u;)(k,l\] . (23) 
k,k' Jo 

To show that Hjy generalizes V, of Eq.(2), it is enough to replace the coupling 
9kk'(w) by g(u) for all values of k and use the identity J2 k \k){k\ = 1. 

In the following, we shall use g(cu) = V(u), for which H IV coincides with 
the standard Friedrichs model interaction V. 

To obtain the solution of the eigenvalue problem 

(H — E)^(E) = , (24) 
where H = Hi + H n + H in + V, we write ^(i?) in its most general form, as 

/oo 
duME,co)\k,oj). (25) 
re re re - 

The action of the Hamiltonian on ^(i?) yields 

(H - E)V(E) = ME)(ck - E)\l) + £ K^mk) 

_ POO 

+ ^2(c k + u -E)(l> ktl (E)\k,l) + ^2 duME,w)V(u)\k,l) 
k k Jo 

poo poo 

+ / dLuV*(Lu)\k,uj) + V / dw^ fc (^,f;)(c fc + w-S)|A;,a;) 

re J ° k J ° 

/oo 
dw^ k (E, uj) /* (to) | Z> + (E) h ik | i, 1) 
re,t i,fe 
/•oo 

+ V/ tfc;^ fc (E)/ ifc ( W )|i,w) = (26) 
i,* ^ 

The above equation can be re-arranged as a linear combination of \k), \k,l), 
and \k,u>) with vanishing coefficients. Thus, 

POO 

ME)(c k - £) +J2hME) + £ / dcuME,u;)ft k (cu) = (27) 
, , Jo 



Vi{E)h kl + (c fc + too ~ E)<p kl (E) + / duME, uj)V(uj) = (28) 
i J o 

ME,u)(ck + u-E) + J2w( E )fki(u)+<f>ki(E)V*(ou) = (29) 

i 

for \k), \k,l), and \k,u), respectively. A new set of equations, where only 
the amplitudes <fi(E) and <pn(E) appear, can be obtained from the above 
system. From the last equation we obtain: 



?pk{E,uj) = c5{c k - u - E) - > — , (30) 

^— ' c fc — u) — E Cfc — u) — E 



where c is an arbitrary constant. 

If we replace (1301) into ( l27j) and (1281 . we obtain an infinite system of 
coupled equations in the amplitudes <fin(E) and fi(E). In the most general 
case, to find the analytic solution of this system may be rather difficult. Note 
that the resonance behavior is attained by the coupling to the bosons. In 
order to obtain solutions for the above system, we need to make appropriate 
choices of the model input represented by the couplings: h kh f k i{uj) and 
V(u). 

After this substitution, we obtain the following result: 



(Cfc - E)5 km - 22 A km(E) 



Vm(E) 



+ iKnk ~ B km (E)) Ki{E) = -cY, flk(E ~ c m ) 



(31) 



and 



J2[hki ~ B kl (E)ME) 



+ ( C fc + u -E- C k (E))(f> kl (E) = -cV(E - c k ) 
where we have introduced the following notations: 



(32) 
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A km (E)= I duo J2 ^ ^ (33) 
" * C\ + uj — hi 



BU E)=I duS^TM (34) 

C m + UJ — Uj 

B tm{ E)=l ^y^J (35) 
' c k + uj — E 







Notice that expressions ( |34ll36l) contain both f(uj) and V(u>), as a result of 
coupling of fermions and bosons in the continuous sector of the model. The 
solution already found in Sections 2.1 and 2.2 is recovered if fkmiu) = h km = 
0, since C k (E),^0 and then Eq.([32]) reduces to Eq.©. 



3 A simplified version of the Model: One 
fermion level. 

In this section, we are going to discuss a simplified version of the Friedrichs 
model including the following ingredients: 

i. ) A single fermionic bound level. This means that there is only one 
term in the sum that appears in Hj in ( 1201) . This term is denoted as 

H F = e\i)(i\, (37) 

where we write \i) to denote the fermion level and thus avoid confusion with 
the boson bound state denoted as |1). 

ii. ) A standard Friedrichs model as described in the previous section. 

iii. ) An interaction between the fermion level and the boson field given 
by the following reduced version of Hjjj: 

H UI = J du[f{u)\i){i,u\ + f*{uj)\i,u){i\}. (38) 
Thus, the total Hamiltonian is given by 
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/■oo 

H — u)q |1)(1| + / da; + e 

Jo 

POO 

+A / [V(u)\u)(l\ + V*(u)\l)(u\]d(v 
Jo 

+ /" du[f(oj)\i)(i,uj\ + r(uj)\i,uj)(i\}. (39) 

The interaction produced by H is depicted in Figure 1. 

Take now the part of the Hamiltonian that corresponds to the standard 
Friedrichs model, given by Hq + XV as in ([1]) and (j2J). It has been shown Q25J) 
that this Hamiltonian can be diagonalized in terms of the Gamow vectors 
and the background part as 

H=z R 1/oX/ol + J z\z)(z\ dn(z) (40) 

In the description given by fj40|) . the boson bound state has been replaced 
by the Gamow vector |/o) and the boson field by the background. We can 
ignore or neglect the background term and then consider a two level system: 
one level is given by the fermion bound state \i) and the other by the Gamow 
vector |/o). The resulting model can be seen as depicted in Figure 2. 

The situation now as visualized in Figure 2 contains two systems in in- 
teraction for which the Hamiltonian can be written aa^| 

H = s\i)(i\ + z R \fo)(fo\+A(z R )[\i,fo)(i\ + \i)(iJo\], (41) 

3 As a matter of fact, this Hamiltonian should be written as 

H = e\i)(i\ ®I + z R I® I/oX/oI + A(z R )[\i, f )(i\ + f \ } , 

where I is the identity operator. In terms of creation and annihilation operators, this 
Hamiltonian can be written as 

H = ec\ Cl + z R G^G + A(z R )[clG^c t + c\G Cl ] , 

where c\ and Ci are the creation and annihilation operators for the fermion with state \i) 
and G^ and G are the creation and annihilation operators for the Gamow state |/o). 
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where the symbol A(z R ) denotes a complex number which depends on the 
position of the resonance (and therefore on the value of the coupling constant 
A) only. 

It is important to remark that in we have neglected the contribution 
of the background integral, (a term of the form f r z\z)(z\ dfi(z)) and note 
that consequently the resulting Hamiltonian is not formally self-adjoint. 

Our goal is to show the existence of resonances, i.e., complex solutions in 
r] of the eigenvalue equation 

Hip = rjip . 

The most general form of the vector ip can be given as: 

<p = a\i,0)+/3\i,f Q ), (42) 

where: i.) The vector \i, fo) is the tensor product \i) <8> |/o). ii.) The vector 
\i,0) is the tensor product \i) ® |0), where |0) denotes a state with zero 
resonances. In the basis given by 0), \i, fo)} the matrix elements of H 
can be readily obtained: 

/ (i,0|#|«,0) (i,0\H\i,f Q ) \ ( e A(z R ) \ 
H= I _ = • (43) 

\ <i,/o|iEr|i,0> (i, f \H\i, f ) J \A(z R ) e + z R J 

In order to obtain the eigenvalues of H, we need to solve the equation 
det(H — E) = and this gives: 

det(H - E) = (e - E)(e + z R - E) - A 2 (z R ) = , (44) 

i.e., 

E 2 - (2e + z R )E + z R e - A 2 (z R ) = , (45) 

which gives 

E = e+ Z f±yz 2 R + A*(z R ). (46) 

Observe that, for any given resonance z R in the boson levels, there exists two 
new resonances in the coupled boson-fermion system. 
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3.1 Calculation of A(zr). 



In the simplified version under our consideration, only the interaction Hamil- 
tonians i/jjj and Hjy remain. For the former, we keep only one term in the 
sum in (1221) so that it has now the form: 



H. 



in 



[f{u)\i,0){i,u\ + r{u)\i,u){i,0\]dto. 



(47) 



Concerning Hyi, it only survives a term of the sum in (I23p and the corre- 
sponding form factor (the g(uj)) is the form factor of the boson-boson field 
coupling V(uj). As it was previously clarified, this coupling is the responsable 
for the creation of the boson resonance given by |/o). 

We now want to determine the numbeiQ A(zr) in the approximation that 
neglects the background. Note that (jUJ) is, save for the background term 
equal to Hj + Hjy + Hjjj, with only one fermion bound term in the fermion 
sum in Hf. Since e\i) (i\ + z R \f ) (/ | + background = H I + H IV , we have that 



Hui = A(z R )[\i,f o )(i,0\ + |z,0)(z,/ |] ■ (48) 
If we multiply ( JiSj) to the right by the ket \i, f ), and take into account that 
(*,0|i,/ ) = and (i,f \i,f ) = (i\i)(fo\fo) = 1 [23 US], we obtain: 



f(u)\i)(i,u\i,f Q )du) 
which obviously gives: 



A(^) 



f{u)(u\f ) duj 



LJQ 



\i) = [A(z R )} \i) , (49) 



/M(w|/o) du. 



(50) 



The value for the bracket (u)\f ) is readily obtained from (fT8"j) . This gives: 



' z R - u + lO 



APV / cko — zttajIzr) V(zr) 

J o zr-uj 



(51) 



4 We write it as a function of the position of the resonance pole, although we equally 
can write it as a function of the value of the coupling constant A which determines the 
pole, i.e., as A(A). 
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Here, PV stands for principal value (although since z R is complex the prin- 
cipal value is here the value of the integral). We are always assuming that 
both form factors, f{oS) and V(u), are analytically continuable and that have 
a certain well defined value at z R . 

Neglecting the integral term in fl5T]) is an approximation similar to ne- 
glecting the background term. From this point of view, we arrive at 

A(z R ) « -m\f{z R )V{z R ) . (52) 
Note that, after (1521) . A (z R ) is complex. 

3.2 The effect of multiple fermion levels. 

We are going to introduce a generalization to the above discussed simplified 
model. Since we are discussing decay in nuclei, it seems natural to add 
multiple levels for the fermion bound states. Now, in the extended Friedrichs 
model we keep the following terms: 

1. - The free Hamiltonian Hj as in (1201) . 

2. - The interaction between fermions and the boson Hjj as in (1211) . 

3. - The interaction between fermions and the boson field Hm given by 
(|22} . 

4. - The interaction boson, boson field given by (J2j. Again, note that this 
is a simplified form of Hjy. 

The strategy is the same as in the previous subsection: we replace the 
boson-boson field by the resonance-background and we neglect the back- 
ground field keeping the Gamow state. Then, we have an interaction between 
the bound state fermions and the boson resonance given by0 

#int = J2 A M^R)[\k,fo)(lA + MXfcJoN ■ (53) 

kl 

5 If c k and Cfe respectively are the creation and annihilation operators for the k — th 
fermion and G^ and G the creation and annihilation of the Gamow, then, this interaction 
Hamilonian can be written in terms of these operators as: 

H iBt = J2 A kiclc k (G^ + G). 

kl 

We can operate either in the language of vectors, as shall do in the main text or in the 
language of creation and annihilation operators as is usual in the second quantization 
formalism. 
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Repeating the procedure of the previous subsection, we identify the in- 
teraction Hamiltonian, excluding AV^, with H int given in (1531) . This gives the 
following identity: 



#int = £[/ifc,|Z,0><fc, 1| + Ki\k, 0| ] 

kl 

poo 

+ Y, [f k i(u>)\l,0)(k,u\ + f* kl (u)\k,u)(l,0\]. (54) 
ki Jo 

Multiplying (1341) to the right by \m, fo), we obtain: 

fki(uj)\l,0)(k 



= J2 A ki\l,0)(k,f \m,f ). (55) 

kl 

Since (l|/o) = 1 and (k\m) = 5k m , (1331) yields to 

POO 

Y\{h m i + f ml (co)(u\fo)du}\l,0) = y] A m ,(zo)M). (56) 
i ^ « 

As the |/, 0) are linearly independent, (1361) gives for all m, I: 

POO 

A m /(^o) = h m i+ / /^(^(wl/o) ^ 

./O 

& h m i-iirf m i(zo)W(z ) , (57) 

where again we have used ffl~8]) for the determination of the value of (oo\fo) 
and have neglected the term principal value as in (I32|) . 

The space of pure states is formed by the vectors of the fornj§ 



It we use the second quantization language, the vector can be written as <p(E) 
Sfc a fc c !|0) + Yl PkCkG^\Q), where |0) represents the state of the vacuum. 
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<p(E) = J2®k\k,0) + J2Pk\k, f ) . (58) 

k 

Then, we want to find the explicit form of the Hamiltonian in terms of the 
basis given by the \k, 0) and the \k, fo) and then obtaining the eigenvalues of 
this Hamiltonian. The presence of complex eigenvalues of this Hamiltonian 
will be equivalent to the presence of resonances for the model. 

In order to show in detail the calculation, we split the Hamiltonian into 
three terms as follows: 



H 1 = ^2e\k,0)(k,0\, (59) 
k 

H 2 = z R \0,f )(oJ \, (60) 

H 3 = J2^ki\k,f )(lJo\ + H BG , (61) 

kl 

where Hbg denotes the background part of the Hamiltonian that we neglect 
and we have written instead of Aki{zo) in order to simplify the notation. 
After simple manipulations, we obtain the following^: 



k' 

(i P \H 2 \ ( p) = J2{P*Mz R , 
w 

(ip\H 3 \<p) = obAv(0) + #M«) , (62) 

k' k' 

where 

7 We arrive to the same results using second quantization notation. In this notation, 
#i = J2k £ k c { c k, H-2 = ZrG^G and H 3 = J2ki^ki{G^ + G)c\ci. We use the anti- 
commutation relation {ck,cj} = Ski and the commutation relation [G, G^] = 1 and note 
that the action of and G^ in the vacuum |0) is given by c^|0) = \k, 0) and G* |0) = |0, fo). 
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Afc'(a) = ^ Afc'fc"a fc » , Afe'(/3) = 2j ^k'k"Pk" • (63) 



fc" fc" 
Using (|62|) and (|43|) . we finally get: 



kk 1 



In order to obtain solutions for the eigenvalue equation Hip = Eip, we 
use the following variational principle: 



(/ {{<p\H\<p)-E(<p\<p)} = (65) 



da* k , 




-{{ip\H\ V )-E{ V \ V )} = Q. (66) 
°Pk> 

This gives a system of linear equations in the undetermined a^, (3^, 
which gives the eigenvector \(p) and E. This system is (we are assuming that 



OLk'^k' + ^2 Akk 'P k ' = Eak ' ( 67 ) 
k' 

Pk> {£k> + zr) + ^2 A kk>a!k = E(3 k > . (68) 

k 

This can be also written as: 



(e v — E) = — ^2 ^kk'Pk , (69) 



A' 



/3fc'(gfc' + - Jg) = - ^ Afcfc'Qfe ■ (70) 



A' 
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Thus, the eigenvalue equation can be written in the following matrix form: 



£ A 

A £ + Z 



with 



a 
b 



E 



(71) 



A 



/ An 

A 2 i 



V A 



Ai 2 
A 22 

kri 1 kn 2 



Ai fcn \ 

A2fc n 

Ai. u I 



£ 



£fc 2 



\ 



\ 





(72) 



£ + Z 



( £ki + zr 




V 







/ Oil \ 

a 2 





£fe 2 + Z R 


/ A \ 

02 








(73) 



(74) 



\ «fcn / \ An / 

where fc n is the total number of fermions considered. 

If we have one fermion state, only, we are in the case studied in the 
previous subsection. If we have two fermions states, the eigenvalue equation 
is a equation of fourth order and therefore solvable. In the simplest case in 
which 



a 
a 



(75) 



a being a complex constant, the eigenvalues of equation (ITT1) fulfill the fol- 
lowing equation: 



det 



/ ei - E a \ 

£ 2 - E a 

a £1 + z R - E 

\ a £ 2 + z R -E ) 



{(ei - E){£ 2 + z R -E)- a 2 }{(£ 2 - E)( £l + z R - E) - a 2 } 



(76) 
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The solutions of (175]) are 



Ei,2 = 7j (z R + £ 1 + e 2 ± V( Z R + £i + ^) 2 - 4(ei£ 2 + e 2 z R - a 2 ) J (77) 

E 3 ,4 = ~ + ei + £2 ± + ii + e 2 ) 2 - 4(eie 2 + - « 2 )) • (78) 
Obviously, these four solutions are complex. 

3.3 Realistic Applications 

In this subsection we shall rephrase some results of the previous section in 
terms of the particle-vibration coupling scheme [I] and, particularly, in the 
context of the particle- vibration coupling in Bergreen's basis [T8] . 

To start with, we shall discuss the structure of the particle-vibration 
coupling, in a resonant basis and by using separable two-body interactions. 
The boson sector of the Hamiltonian is usually treated in terms of the TDA 
(Tamm-Damcoff) approximation, leading to vertex functions of the type [H] 

where V is the matrix element of the two-body interaction, v and 77 denote 
particle-hole configurations and u is the energy of the one phonon state. If 
in the configuration 77 (or v) a single-particle resonant state participates, the 
energy of the one phonon state, 00, becomes complex and the imaginary part 
of this energy will be proportional to the imaginary part of the energy of the 
particle-hole configuration. Tipically, in a large single particle basis, with few 
resonant states, the imaginary part of the energy of a high-lying one-phonon 
state is of the order of 1 keV or smaller [18]. 

In the context of the results of the previous subsection, there is a direct 
correspondence between Eq. 51 and A(a>, 77). Note that, while A(zr) has been 
obtained by working with the fermion-boson interaction firstly, the result 
A(u;,77) has been obtained (see [IB]) by diagonalizing the boson sector and 
then by using it in the transformation of the particle-vibration coupling. 

The just mentioned analogy suggest that both procedures may be equiv- 
alent. These procedures are: 
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i. ) The treatment of particle- vibration couplings in resonant (Bergreen) 
basis. It consists in including fermion resonant states in single-particle basis, 
solving for TDA or RPA equations and getting complex phonon energies, 
extracting vertex functions and performing a particle- vibration coupling cal- 
culation between effective fermion and boson fields. 

ii. ) The extended Friedrichs model: It consists in solving the Friedrichs 
model equations proposed in this work, isolating a boson resonant state, 
ignoring the background and performing a particle- vibration coupling scheme 
by using bound fermions and the resonant boson. 

In the case of the simplified versions of section 3.2, the solution E4, in 
(78), closely resembles some of the results of [T5]. For the case of one particle 
above the Fermi level £4 yields a physical value 

E A tte x -K\z R ) (80) 

where we have set a = A(zr). Thus, the coupling with the resonant boson 
state leads to a imaginary energy of the order of the imaginary part of the 
quantity A(zr). 

4 Conclusions 

In this work we have presented an extended version of the Friedrichs model. 
It consists of couplings between fermion and bosons, in addition to the cou- 
plings between the discrete and continuous boson excitations. The extended 
model is exactly solvable and the structure of the solutions shows a suitable 
mechanism for the appearance of complex energies in the fermion spectrum. 
In order to show the possible relevance of the present results for nuclear 
structure physics we have performed a comparison with the results obtained 
in realistic coupling schemes. The coincidence of the formal structure of the 
solutions, for fermion states, in both approaches, suggests a direct conec- 
tion between them. The results allow for the interpretation of the imaginary 
part of the energy of the fermion states obtained by coupling fermion bound 
and resonant states to collective nuclear particle-hole excitations. Further 
studies, concerning the calculation of two-particle correlation functions and 
particle-phonon correlation functions in resonant basis are in progress. 
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